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1 Introduction 



Although the Standard Model of electroweak interactions has worked pretty well so far, it 
has a number of shortcomings. For example, it is difficult to explain why the Higgs mass is 
much smaller than the Planck scale (known as the "hierarchy problem"), a description of 
gravity is lacking and there is no good candidate for dark matter. Moreover, given the fact 
that the Higgs mass is either relatively light or rather heavy, i.e. 122 GeV < Mjj < 127 GeV 
or Mjj > 600 GeV at 95% confidence level with a possible Higgs discovery at a mass 
near 125 GeV [1-5], it is difficult to guarantee the stability of the Standard Model Higgs 
mechanism up to the Planck scale [6]. This might hint at the possibility that there is a 
scale of new physics between the electroweak scale and the Planck scale. 

If signs of such new physics are observed at the upcoming runs of the Large Hadron 
Collider (LHC), the main theoretical challenge will be to unravel the underlying theory. 
Since the underlying physics can reside at energies that largely exceed the reach of the 
LHC, this is not going to be a simple task. The standard strategies for addressing this 
issue make use of renormalisation group techniques, linking the "low-energy" physics ob- 
served at present-day experiments to high-energy parameters at the energy scale where 
the underlying theory is formulated. The most widely studied approach is the top-down 
method, where one starts by choosing a specific new-physics model. Subsequently, the 
high-scale model parameters are evolved down to the collider scale and predictions can be 
made about the way the model will manifest itself phenomenologically in ongoing experi- 
ments, allowing a confrontation between theory and experiment. This approach is great for 
identifying where one should look for signs of new physics, but it is not really well-suited 
for deriving conclusive statements about the underlying model. Alternatively a bottom-up 
method can be used, where one is guided by experimental data. One starts by adopt- 
ing a rather general phenomenological framework for describing the physics beyond the 
Standard Model, such as a super symmetric extension of the Standard Model. Within this 
phenomenological context, the data are converted into running parameters at the collider 
scale. Subsequently, the running parameters are evolved up towards the scale where the 
underlying physics is presumably residing, allowing the high-scale parameters to be con- 
fronted with specific predictions from new-physics models. This method is better-suited 
for getting information on the underlying model. However, since the renormalisation group 
evolution is a numerical procedure, it can lead to an increase of the error from one scale 
to another. Moreover, the fact that we do not know the scale of new physics can lead to 
misinterpretation of the data. A more detailed discussion of the advantages and drawbacks 
of the top-down and bottom-up methods will be given in section 2. 

Recently, an elegant third approach to probing high-scale physics has gained some 
interest [7-10]. This approach is based on the same philosophy as the bottom-up method, 
but instead of using the full set of running parameters it makes clever use of so-called 
renormalisation group invariants. These are combinations of running parameters chosen in 
such a way that they do not evolve with energy. If we measure their values at the collider 
scale, we will immediately know their values at the threshold of new physics. This fact 
allows one to probe physics at high energy scales without having to evolve all parameters. 
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Subsequently, the renormalisation group invariants can be combined into sum rules that 
test the underlying physics. Up to now, sum rules have been constructed for testing specific 
models. We advocate to employ sum rules in a more model-independent way, by using them 
as fast diagnostic tools to test generic properties that are common to new-physics models, 
such as unification and universality properties. If a certain property is realised in Nature, 
all corresponding sum rules must be satisfied. So, the main strength of invariant sum rules 
is their falsifying power. 

In order to give an idea of how these renormalisation-group- invariant sum rules work in 
practice, we work them out in detail for the study of supersymmetry-breaking mechanisms 
in the context of a phenomenological version of the Minimal Supersymmetric Standard 
Model. This model is one of the prime beyond-the-Standard-Model frameworks to be 
tested at the LHC, since it offers solutions to several of the problems that plague the 
Standard Model. 

The structure of the paper is as follows. In section 2 we describe the concept of 
effective field theories and discuss renormalisation group techniques. In section 3 we give 
the salient details of the Minimal Supersymmetric Standard Model and a few popular 
supersymmetry-breaking mechanisms. Subsequently, the one-loop renormalisation group 
invariants are listed for a phenomenological version of the model. In section 4 we give a 
detailed discussion of the model-independent and model-specific sum rules that can be used 
for studying the supersymmetry-breaking mechanisms. We will conclude in section 5. 

2 Effective Lagrangians and renormalisation group equations 

In a renormalisable quantum field theory, all divergences can be dealt with by redefining the 
masses and couplings, which become scale dependent. On the other hand, if the Lagrangian 
contains operators with dimension greater than four, an infinite number of counterterms 
is needed to remove all divergences and all predictive power is lost. Hence, the prevalent 
view on renormalisation used to be that a sensible theory describing Nature could only 
contain renormalisable interactions. This view has changed since Wilson's work on the 
Renormalisation Group (RG) [12, 13]. The use of a momentum cutoff A is now considered 
to have a physical meaning: A is the scale at which new physics becomes relevant. For 
processes at energies greater than A, the theory is not valid anymore and should be replaced 
by a more fundamental theory. This is the motivation for using effective field theories 
(EFTs). 

2.1 Effective field theories 

Physics problems usually involve widely separated energy scales, which allows us to study 
low-energy dynamics without needing to know the details of the high-energy interactions. 
Effective field theories (see e.g. [14]) are the theoretical tool to describe low-energy physics, 
where 'low' means low with respect to some energy scale A. An EFT only takes into 
account states with mass m« A; heavier excitations with m>A are integrated out from 
the action. The information about the heavy states is then contained in the couplings 
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of the low-energy theory: we get non-renormalisable interactions among the light states, 
organised as an expansion in powers of energy/A. 

An effective field theory is characterised by some effective Lagrangian: 

£ = J>0;, (2.1) 

i 

where the Oi are operators constructed from the light fields and the q are couplings contain- 
ing information on any heavy degrees of freedom. Since the Lagrangian has dimension 4, 
dimensional analysis yields: 

Pi] = di c l ^ 7 J- I , (2.2) 

where A is some characteristic heavy scale of the system. At low energies, the behaviour 
of these operators is determined by their dimension: 

• Operators with di < 4 are called relevant, since they give rise to effects that become 
large at low energies. 

• Operators with di > 4 are called irrelevant: at energy scales E their effects are 
suppressed by powers of E/A, making them small at low energies. These are non- 
renormalisable operators that contain information about the underlying dynamics at 
higher scales. 

• Operators with di = 4 are called marginal, because they are equally important at all 
energy scales. 

This explains why we are able to include non-renormalisable operators in an EFT without 
spoiling its predictive power: at low energies E, their effects can be either neglected or 
incorporated as perturbations in powers of E/A. At high energies, it is more appropriate 
to use a different EFT. Thus at sufficiently low energies, an EFT automatically contains 
only renormalisable operators. 

2.2 Matching 

Suppose we have two EFTs: one that includes a heavy particle and one where its effects 
are included in the form of higher-dimensional operators, suppressed by inverse powers of 
the heavy particle mass M. Since physics around the mass scale M should not depend on 
our choice of theory, both EFTs should yield the same physical predictions. Hence they are 
related by the matching condition: at the threshold [i = M, the two EFTs should give rise 
to the same S'-matrix elements for light-particle scattering. This leads to relations between 
the parameters of the high-energy EFT (the one we use above threshold) and those of the 
low-energy EFT (the one we use below threshold). In other words, the matching conditions 
encode the effects of the heavy field into the low-energy EFT parameters. 

As an example, consider the beta decay of a neutron. In the Standard Model, this 
decay is mediated by a W boson with mass M\y. At energies fi <C Myy, there is not enough 
energy available to produce a physical W boson. Hence we might as well integrate out the 
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W field from the action. In the resulting low-energy EFT, the W exchange reduces to a 
non-renormalisable contact interaction. Matching the two EFTs at jx = M\y yields the 
formula for the Fermi coupling constant 

where g2 is the weak coupling constant. Note that although the W field is not included in 
the low-energy EFT, its 'fingerprints' (namely its coupling constant 52 and mass Mw) are 
still present in the low-energy coupling Gf- Also note that the irrelevant operator corre- 
sponding to the contact interaction is indeed suppressed by powers of Mvp, as mentioned 
in section 2.1. 

2.3 How to look for clues about new physics 

In the process of renormalising a theory, we redefine the masses and couplings by having 
them depend on a reference scale [i. This //-dependence is given by the renormalisation 
group equations, which can be calculated from the condition that anything observable 
should be independent of //. These equations depend on the loop order of the calculations 
and on the particle content of the theory. 

In the framework of effective field theories, physics is described by a chain of EFTs. The 
RG equations and matching conditions for these EFTs allow us to evolve a theory between 
different energy scales. Suppose that we have an EFT describing physics at some (high) 
energy scale \i. If we want to describe physics at a lower energy scale, we have to evolve 
down the running parameters using the RG equations of this EFT. Whenever we reach 
a particle threshold, we must switch to an EFT that does not contain the corresponding 
field. The appropriate matching conditions at the particle threshold yield the parameters 
of the low-energy EFT. From there we can continue to evolve these parameters, now using 
the RG equations of the low-energy EFT. 

From this point of view, the Standard Model is only a low-energy EFT of Nature. The 
shortcomings of the Standard Model hint at the existence of a more fundamental theory. 
Even if that more fundamental EFT is appropriate only at energies beyond experimental 
access, the idea of a chain of EFTs certainly helps us study that more fundamental theory 
(figure 1). We could measure the running parameters at a low scale jx and then evolve them 
upwards. At the threshold of the more fundamental theory, the matching conditions act 
as boundary conditions for the renormalisation group. Hence, by comparing our evolved 
masses and couplings with the predicted matching conditions, we can get information on 
the high-energy theory. 

2.4 How to probe the high scale 
2.4.1 Top-down method 

The literature offers various approaches to using the renormalisation group to extract 
information about high-scale matching conditions. In the context of supersymmetry, the 
most widely studied one is the top-down method. It is called this way because a top-down 
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Figure 1. Scheme for studying physics at scales beyond experimental access. The running couplings 
are measured at a scale where the low-energy EFT is applicable. Using the RG equations, they 
are evolved towards the threshold where new fields presumably enter the theory. Then they can be 
compared with the matching conditions predicted by the more fundamental theory. 



study is started from the high scale and the theory is evolved down to the collider scale 
(say O (1 TeV)). One starts by choosing a specific high-scale new-physics model with a 
limited number of free parameters and proceeds as follows: 

• Pick a point in the high-scale parameter space of the model and translate this into 
values of the running parameters at the new-physics threshold (masses, gauge and 
Yukawa couplings, etc.). 

• Evolve the running parameters down to the collider scale using the renormalisation 
group equations. 

• Using the resulting parameter values, calculate the relevant branching ratios and 
cross sections. 

• Compare the results to experimental data and extract constraints on the high-scale 
parameter space of the model. 

The top-down method is suitable for making general phenomenological predictions. For 
example, it is used to find collider signatures that are characteristic for super symmetry. 
However, for the purpose of testing a new-physics model, this method has some serious 
limitations: 

• With the top-down method one can only determine the regions in the high-scale 
parameter space of the model that are consistent with the data. If only small portions 
of the total parameter space seem phenomenologically viable, one might conclude that 
the model is neither likely to be correct nor natural. However, it seems unlikely that 
we can strictly exclude a model this way. 

• In this approach the confrontation between theory and experiment is performed on 
the basis of individual models. Given the fact that there are many high-scale models 



-6- 



on the market and most of them can also be tweaked in one way or another, it seems 
unlikely that we can find the right model this way. 

• Scanning the entire high-scale parameter space of even a single model is very time- 
consuming. To scan it properly, one ought to use a reasonably fine grid and check each 
point separately. But the parameter space is usually too big to perform a full detector 
simulation for each point. For general predictions of supersymmetric phenomenology, 
for example, one usually resorts to using a limited set of benchmark points (see e.g. 
[15]), because many points in parameter space have a very similar phenomenology. 
However, for the purpose of excluding a certain model, this is no satisfying solution. 

• Fitting the numerical predictions to the experimental data becomes much more diffi- 
cult as the number of high-scale model parameters is increased. Therefore one always 
limits oneself to a model with few parameters. But there is no reason to think that 
Nature would restrict itself to only a few parameters in the EFT beyond the next 
threshold. 

2.4.2 Bottom-up method 

The bottom-up method is an alternative to the top-down method. It works by evolving the 
theory upwards from the collider scale to the new-physics threshold. A bottom-up analysis 
consists of the following steps: 

• Convert experimental data into the running parameters at the collider scale. 

• Using the renormalisation group equations pertaining to the low-energy EFT (e.g. the 
Standard Model above the electroweak scale, or the Minimal Supersymmetric Stan- 
dard Model above the supersymmetric mass scale), evolve these running parameters 
towards the scale where new physics presumably comes into play. 

• Analyse the structure of the high-scale parameters: do they fit the matching condi- 
tions predicted by any new-physics model? 

This method seems more suitable for excluding new-physics models than the top-down 
method. Also, there is no practical need to only consider high-scale new-physics models 
with a limited number of free parameters. Furthermore, no time-consuming scanning of 
that parameter space is involved. However, the bottom-up method presents challenges of 
its own: 

• The running parameters at the collider scale will come with experimental errors. To 
determine the uncertainty in these parameters at a higher scale, we also have to 
evolve the error bars. These may become larger while numerically performing the 
renormalisation group evolution, which could make it difficult to tell for example 
whether certain parameters unify or not. 

• We do not know the value of the high scale that should be taken as the new-physics 
threshold; this scale has to be guessed. In practice, one might evolve the running pa- 
rameters until some of them unify and take the corresponding scale as the new-physics 
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threshold. But a unification scale does not necessarily correspond to a threshold. 
Also, there may be an intermediate new-physics threshold even though no unifica- 
tion occurs there. In both cases, we would extract incorrect boundary conditions for 
matching with the underlying high-energy EFT. 

• Because the RG equations are coupled, all running parameters must be known. Hence 
if we fail to measure one mass or coupling, the bottom-up method cannot be used 
except for subsets of parameters whose RG equations contain only the parameters 
from that subset. 

2.4.3 Renormalisation group invariants 

Recently, a third approach to probing the high scale has gained some interest [7-10]. 
This approach makes clever use of renormalisation group invariants (RGIs). These are 
combinations of running parameters chosen in such a way that they are independent of the 
renormalisation scale [i. A well-known example of an RGI is the following combination of 
the strong coupling #3, the weak coupling 52 as well as the scaled hypercharge coupling 

91 =g'V^- 

Ig = (h ~ h)gf + (63 - h)g^ 2 + (&! - 6 2 )g 3 - 2 , (2.4) 

where the gauge couplings satisfy the following renormalisation group equations at one 
loop: 

16tt 2 ^ = V^ (a = 1,2,3). (2.5) 

Here t = In (fi/fio), where is a reference scale that makes the argument of the logarithm 
dimensionless; its value is arbitrary since it drops out of the RG equations. The RG 
coefficients b a are constants depending on the particle content of the model. For the 
Standard Model they are b a = (fjj,— if) - 7) for a = 1,2,3, whereas for the Minimal 
Super symmetric Standard Model (MSSM) they are b a = (P,l,— 3). It is easily checked 
that dlg/dt = 0. Note that I g is not exactly RG invariant, since we used the one-loop RG 
equations to construct it. We will come back to this issue in section 4. 

A crucial property of RGIs is that if we measure their values at the collider scale, we 
will immediately know their values at the threshold of new physics. This fact allows us to 
probe physics at high energy scales without having to evolve all parameters. For example, 
I g can be used to test whether the gauge couplings unify. Note that if the gauge couplings 
have a universal value g v at some energy scale, then I g will vanish. Since I g is an RGI, 
it will consequently vanish at every scale where the renormalisation group equations (2.5) 
are valid. Hence, if we measure the gauge couplings at one scale, we can perform a quick 
diagnostic check to test whether gauge-coupling unification occurs within the context of a 
specific EFT. Also note that this consistency check is independent of the value of the scale 
where the gauge couplings unify. 

To illustrate this, we perform this check explicitly for the Standard Model, using the 
measured couplings at fj, = Mz in the MS scheme. The gauge couplings can be obtained 
from the measured quantities a^ l (M z ) = 8.45 ± 0.05, a~ l (M z ) = 127.916 ± 0.015 and 

1 This occurs for example in a supersymmetry-breaking model called Mirage Mediation, see section 3. 
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Figure 2. Two- loop renormalisation group evolution of the inverse gauge couplings squared 
a a 1 {P L ) — 47T3a 2 (A t ) f° r a = 1,2,3 in the Standard Model (dashed lines) and the Minimal Su- 
persymmetric Standard Model (solid lines). For the latter, the supersymmetric particle thresholds 
are varied between 500 GeV and 1.5 TeV, and a^{Mz) is varied between 0.117 and 0.121. Figure 
taken from [17]. 

sin 2 6 W (M Z ) = 0.23131 ± 0.00007 [16]. Here a" 1 = 4irg~ 2 , a is the fine structure constant 
and 9w is the weak mixing angle. The latter two are related to ct2 and aq by the relations 

a 2 _1 (M z ) = sin 2 9 w (M z )ar 1 {M z ) = 29.588 ± 0.010, 

a7 l (M z ) = - cos 2 e w (M z )a~ 1 (M z ) = 58.997 ± 0.009. (2.6) 
5 

Using these values we find J| M = -3.252 ± 0.030, which lies many standard deviations 
from 0. Hence we find no compatibility with gauge-coupling unification within the Standard 
Model; this can be confirmed using the bottom- up method (figure 2). 

Similarly, we can make an estimate for I g in the MSSM by taking the above values of 
a- l {M z ). 2 We find I^ssm = _o.059±0.024, which is close to zero. Hence the MSSM might 
allow for gauge-coupling unification, depending on the actual values of the supersymmetric 
particle thresholds. This can be confirmed using the bottom-up method (figure 2): the 
gauge couplings are observed to unify, even when we use the two- loop RG equations. 

To summarise, RGIs provide a fast diagnostic tool for probing matching conditions at 
high energy scales. They circumvent the need to evolve the running parameters numerically; 
we do not even need to know exactly at which energy scale new physics arises. 

2 In order to determine the actual value of /MSSM wg OU gj 1 ^ ^ Q use ^ ne values of a„' at the scale where 
the MSSM becomes valid, i.e. the highest supersymmetric particle threshold. This scale will be somewhat 
higher than Mz- As can be seen using the bottom- up method (figure 2), this barely changes the prediction 
of gauge-coupling unification: only the value of the unification scale might change. Hence we should get a 
good estimate of /^ 1SSM using the values of a~ 1 (Mz)- 
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3 Using RGIs and sum rules to study supersymmetry breaking 

In order to give an idea of how the RGI method works in practice, we work it out explicitly 
for models of supersymmetry breaking. To this end, we first give a short description of the 
salient details of the minimal supersymmetric extension of the Standard Model. This is 
followed by a discussion of the necessity to break supersymmetry and the ways to achieve 
this. For an extensive introduction to the subject the interested reader is referred to 
Refs. [17, 18]. 

3.1 Supersymmetry 

Although the Standard Model has worked pretty well so far, it has a number of short- 
comings which we mentioned in the introduction. Supersymmetry offers possible solutions 
to these problems [17, 18]. It is a symmetry between fermions and bosons that would 
complete the list of possible spacetime symmetries [19, 20]. The single-particle states of 
a supersymmetric theory fall into irreducible representations of the supersymmetry al- 
gebra, called supermultiplets. These contain an equal number of bosonic and fermionic 
degrees of freedom. The bosons and fermions in a supermultiplet are called superpartners 
of each other. A minimal supersymmetric extension of the Standard Model contains two 
types of supermultiplets: chiral (matter) supermultiplets, which consist of a two-component 
Weyl spinor and a complex scalar field, and vector (gauge) supermultiplets, which consist 
of a spin-1 gauge-boson field and a spin-1/2 Majorana spinor, called the gaugino field. 
The supersymmetry generators Q, commute with the mass-squared operator P 2 and 
the generators of gauge transformations, so superpartners have the same mass and gauge 
quantum numbers. In view of the quantum-number structure of the Standard Model, this 
implies that a supersymmetric extension of the Standard Model introduces at least one 
new supersymmetric particle (or sparticle for short) for each Standard Model particle. 

3.2 The Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric Standard Model (MSSM) is defined to be a supersymmet- 
ric extension of the Standard Model with minimal particle content and a new conserved 
quantum number called i?-parity [21]. The MSSM particle content and nomenclature is 
listed in tables 1-2. The only new parameter with respect to the Standard Model is a 
supersymmetry-preserving Higgs-mixing parameter \i. 

3.3 Constraints on broken supersymmetry 

If supersymmetry were an exact symmetry of Nature, each sparticle would have the same 
mass as its Standard Model partner and we would have discovered them already. Hence, 
if supersymmetry is a symmetry of Nature, it must be broken somehow. The requirement 
that broken supersymmetry should still solve the problems of the Standard Model puts 
constraints on the possible terms of a supersymmetry-breaking Lagrangian. 

In order to maintain the solution to the hierarchy problem, we must consider soft su- 
persymmetry breaking. This means that we only consider supersymmetry-violating terms 
with masses and couplings of positive dimension. By using only relevant operators (see 
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Symbol 


Spin 


Spin 1/2 


representation 
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Q = (u L , d L ) 




(3,2,i) 


squarks & quarks 
(3 generations) 
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■ ' 5fc 
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(3,1,-|) 
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a R 


d C R 


(3,1,|) 


sleptons & leptons 


L 


L = (v, e L ) 


(y, e L ) 


(1,2,-^) 


(3 generations) 


e 


e R 




(1,1,1) 




H u 






(1,2, 1) 


Higgs & higgsinos 


H d 


(H° d , H-) 


, ^) 


(1,2,-|) 



Table 1. Chiral supermultiplet content of the Minimal Supersymmetric Standard Model and the 
corresponding representations of the gauge group SU(3)c X SU(2)l X t/(l)y, with superpartners 
indicated by a tilde. Note that we need an additional Higgs doublet compared to the Standard 
Model and that right-handed modes are charge conjugated in order to bring them into left-handed 
form. 



Names 


Spin 1/2 


Spin 1 


Gauge- group 
representation 


gluino & gluon 


9 


9 


(8,1,0) 


winos & W bosons 


w 1 w 2 w 3 


w 1 W 2 W 3 


(1,3,0) 


bino & B boson 


B° 


B° 


(1,1,0) 



Table 2. Gauge supermultiplet content of the Minimal Supersymmetric Standard Model and the 
corresponding representations of the gauge group SU(3)c x SU(2)l x U(1)y- 

section 2.1) to break supersymmetry, we guarantee that the high-scale physics responsi- 
ble for supersymmetry breaking decouples at low energies. It also guarantees that non- 
super symmetric corrections to the Higgs mass vanish in the limit m so ft — > 0, where m so ft 
is the largest mass scale associated with the soft parameters. 

The most general soft-supersymmetry-breaking Lagrangian in the MSSM, compatible 
with gauge invariance and -R-parity conservation, contains complex gaugino masses Mi, 
M2, M%; trilinear couplings a u , a d and a e , which are complex 3x3 matrices in family 
space similar to the Yukawa couplings; sfermion mass terms mq,mQ,m|, m^mf, which 
are Hermitian 3x3 mass matrices in family space; real Higgs masses mn 2 u u and m 2 ^ and a 
complex supersymmetry-breaking Higgs mixing parameter b. 
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Soft super symmetry breaking introduces 97 new masses, mixing angles and phases 
[22]. These may lead to predictions for flavour-changing neutral current (FCNC) and CP- 
violating processes that violate experimental bounds. In order to suppress these processes, 
we additionally assume the following relations between the parameters: 

• The soft sfermion masses are flavour diagonal and the first- and second-generation 
masses are degenerate. 

• There are no sources of CP-violation in the soft-supersymmetry-breaking sector be- 
yond those induced by the Yukawa couplings. 

We also neglect the first- and second-generation Yukawa and trilinear couplings, because 
they give very small contributions to the evolution of the soft-supersymmetry-breaking 
parameters. These contributions are smaller than the two- loop corrections associated with 
the gauge couplings and third-generation Yukawa couplings. 

These assumptions resemble the ones that form the basis for the so-called phenomeno- 
logical MSSM (pMSSM) [23]. Note, however, that we opt to work with the soft-supersym- 
metry-breaking parameters in our approach rather than the mass eigenstates that are used 
in the pMSSM. Under these assumptions, we are left with the following parameters: 

• Twelve real soft scalar masses, which we denote as m~ , m~ , m~ , m~ , m~ , m~ , 

Qi' Qs' «i' «s' di' d 3 

m~ , m~ , mk , mA , m|r , m 2 u , in accordance with the notation in table 1. The 

Li' L a ei> e 3 > ««' «d' 

subscripts 1 and 3 refer to the first and the third generation respectively. 

• Three real gauge couplings g±, 52 , 93- 

• Three real gaugino masses Mi, M%, M3. 

• Three real third-generation Yukawa couplings yt, yb, Ut- 

• Three real soft trilinear couplings At, Aj,, A T defined by 

di = Aiyi (i = t,b,T, no summation). (3-1) 

• Two real Higgs mixing parameters fi, b. 

One might even go one step further and assume full universality, i.e. take all mass 
matrices proportional to the unit matrix. Such apparently arbitrary relations between 
the soft parameters could make sense from the effective-field-theoretical point of view. If 
supersymmetry is exact in a more fundamental EFT than the MSSM, but is broken sponta- 
neously at some high energy scale, then the terms in £ S oft may arise as effective interactions. 
In that case, the universality relations could arise as matching conditions at the thresh- 
old where we switch from the more fundamental theory to the MSSM. Strictly speaking, 
flavour-universality is lost once the parameters are evolved down to the electroweak scale, 
but the numerical impact of this RG evolution is small [24]. Thus the desire for a theory 
that naturally explains supersymmetry breaking forces us to consider spontaneously broken 
supersymmetry. 
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3.4 Breaking supersymmetry 

For a spontaneous breakdown of supersymmetry we need a Lagrangian that preserves 
supersymmetry but a vacuum state that breaks it. During the construction of a super- 
symmetric theory, one has to introduce auxiliary fields for each supermultiplet in order to 
make the supersymmetry algebra close off-shell. These are scalar fields that turn out to be 
suitable for breaking supersymmetry (see e.g. section 7 of [17]): if some of them acquire a 
non-zero vacuum expectation value (VEV), supersymmetry is broken. 

It turns out to be difficult to make this happen using only renormalisable interactions 
at tree level. Therefore, the MSSM soft terms are expected to arise radiatively. In radia- 
tive supersymmetry-breaking models, supersymmetry is broken in a hidden sector, which 
contains fields that have no direct couplings to the MSSM fields. The latter are said to 
be in the visible sector. The two sectors only interact indirectly; the interactions between 
them are responsible for mediating the supersymmetry breakdown from the hidden sector 
to the MSSM. If the mediating interactions are flavour blind, then the soft terms of the 
MSSM will automatically satisfy universality conditions. We will discuss several proposals 
for breaking mechanisms in the sections 3.5-3.8. 

3.5 Supergravity 

Supergravity (SUGRA) [25] is the theory that results from imposing local supersymmetry 
invariance. Recall that once we promote a global gauge symmetry (with bosonic gener- 
ators, satisfying commutation relations) to a local one, we have to introduce a bosonic 
field with predetermined gauge-transformation properties. Similarly, by promoting su- 
persymmetry (which has fermionic generators, satisfying anizcommutation relations) to 
a local symmetry, we have to introduce a fermionic field ^ with spin-3/2. This is the 
gravitino, the superpartner of the spin-2 graviton. The resulting SUGRA Lagrangian is 
non-renormalisable; there is as yet no renormalisable quantum field theory of gravity. How- 
ever, the non-renormalisable operators are suppressed by inverse powers of the Planck mass 
Mpi = 0(1O 19 GeV), so that their effects at low energies are small (see section 2.1). 

The spontaneous breakdown of supersymmetry occurs in a hidden sector where the 
auxiliary component of some superfield gets a VEV. According to Goldstone's theorem, 
spontaneously breaking a global symmetry yields a massless particle with the same quan- 
tum numbers as the broken symmetry generator. Since the broken generator Q is fermionic, 
the massless particle is a massless neutral Weyl fermion, called the goldstino. The gold- 
stino then becomes the longitudinal component of the gravitino, which becomes massive. 3 
It turns out that when we consider the effects of the supersymmetry-breaking VEV, the 
gravitino mass sets the scale of all the soft terms. Moreover, the scalar masses are 

universal at the scale where supersymmetry becomes broken. 

Minimal supergravity. The most widely used model of supersymmetry breaking is min- 
imal supergravity (mSUGRA) [26, 27]. Despite the name, mSUGRA is not a supergravity 
model, but rather the low-energy EFT resulting from a minimal locally super symmetric 

3 Because of the similarities with the Higgs mechanism, where the electroweak gauge bosons 'eat' the 
Goldstone bosons and become massive, this mechanism is called the super-Higgs mechanism. 
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model. In the underlying model, one uses the simplest possible Ansatz for the scalar po- 
tentials. This leads to universal soft supersymmetry-breaking parameters in the scalar 
sector. Gauge-coupling unification in the MSSM suggests an additional simple Ansatz for 
the gauge kinetic function, which leads to universal gaugino masses. The model is then 
described by four parameters and a sign. This involves three parameters at the GUT scale 
-^GUT = 2 • 10 16 GeV: a universal scalar mass mo, a universal gaugino mass Mi/ 2 and a 
universal proportionality factor Aq between the trilinear couplings and the corresponding 
Yukawa couplings. In addition it involves the ratio tan f3 of the two non-zero Higgs VEVs 
and the sign of the supersymmetric parameter /i, which are both usually determined at the 
weak scale. At the GUT scale, the soft terms relevant to our study are therefore given by 

mf = uiq, (3.2a) 
M a = M 1/2 (a = 1,2,3), (3.2b) 

where mf are the scalar squared masses. From the supergravity point of view, the pa- 
rameters mo, M1/2, A) depend on the hidden-sector fields and are all proportional to m 3/ / 2 
(for example, one has the relation mo = m 3/ / 2 ). However, from the perspective of the low- 
energy EFT that we call mSUGRA, they are simply regarded as model parameters. The 
MSSM is assumed to be valid up to the GUT scale, where the relations (3.2) serve as RG 
boundary conditions. In addition, the soft Higgs mixing term B = b/fi has the GUT-scale 
value Bq = Aq — m 3 / 2 . 

As an aside, there is a model similar to mSUGRA: it is called constrained MSSM 
(CMSSM, see e.g. [28]). It has the same boundary conditions as mSUGRA and these two 
models are often confused in the literature. However, mSUGRA arises from a supergravity 
model whereas the CMSSM does not: the boundary conditions (3.2) are simply postulated. 
Also, in the CMSSM there is no relation between the model parameters and m 3 / 2 , and the 
relation Bq = Aq — m 3/ / 2 does not hold either. 

Supergravity models are attractive since they provide a natural framework for super- 
symmetry breaking: a locally supersymmetric Lagrangian automatically contains terms 
that can mediate supersymmetry breaking. In addition, mSUGRA has great predictive 
power since it has only four free parameters plus a sign. However, these models must 
necessarily appeal to Planck-scale physics, which is still poorly understood. Furthermore, 
though gravity is flavour blind, the supergravity invariance of the Lagrangian cannot pre- 
vent the occurrence of (Planck-scale suppressed) flavour-mixing operators that correspond 
to tree-level interactions between hidden-sector fields and visible-sector fields. In order 
to suppress sparticle-induced FCNC processes, one must resort to additional generation 
symmetries. 

3.6 Anomaly-mediated supersymmetry breaking 

In some models of supergravity, the visible and hidden sectors are physically separated by 
extra dimensions [29, 30]. In these 'braneworld' scenarios, often inspired by string theory, 
our four-dimensional world is embedded in a higher-dimensional bulk that has additional 
spatial dimensions, which are curled up. 
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The general idea is that the MSSM fields and the hidden-sector fields are confined 
to parallel, distinct three-branes (space-like hypersurfaces) , separated by a distance r. 
Only the gravity supermultiplet (and possibly new heavy fields) resides in the bulk. In 
this scenario every flavour-violating term that plagues supergravity, caused by tree-level 
couplings with a bulk field of mass m, is suppressed by a factor e~ mr . Provided that r is 
large enough, the flavour-violating effects are exponentially suppressed without requiring 
any fine-tuning. This class of models is called Anomaly-Mediated Supersymmetry Breaking 
(AMSB), because the size of the soft supersymmetry-breaking terms is determined by the 
loop-induced superconformal (Weyl) anomaly [31]. Local superconformal invariance is a 
rescaling symmetry that is violated at the quantum level. 

Anomaly-mediated terms are always present in supergravity, but they are loop-sup- 
pressed with respect to the gravitino mass and hence result in subleading-order contribu- 
tions to the soft masses. AMSB is the scenario where there are no supergravity contribu- 
tions at tree level, so that the anomaly-mediated terms become the dominant ones. At the 
scale M AMSB where supersymmetry breaking occurs, the soft terms relevant to our study 
have the following values (using the usual pMSSM assumptions): 

M a = 9 2 a (M AMSB )m 3/2 (a = 1,2,3), (3.3a) 

m \ = \ii m \li- ( 3 - 3b ) 

Here, mf are again the scalar squared masses with 73 being the corresponding anomalous 
dimensions and b a = 1, — 3) for a = 1, 2, 3. The derivatives 7, = d 7 j/di are explicitly 
given by 
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where we have defined the following quantities for convenience: 

Bt = 6 1 £/* 1 2 + \Vb\ 2 ~ yffl - 3 d ~ j^9i, (3.5a) 

B b = Q\yt\ 2 + \yt\ 2 + \y T \ 2 - ys 3 2 - 3g\ - ^g 2 , (3.5b) 

B T = 4\y T \ 2 + 3|y b | 2 - 3g\ - \gl (3.5c) 

5 
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Figure 3. Contribution to the soft squared mass of the up squark in models of Gauge-Mediated 
Supersymmetry Breaking (GMSB). The auxiliary component of the spurion field X obtains a 
supersymmetry-breaking VEV (Fx)- The up squark only couples indirectly to this VEV: the scalar 
component cj>M of a messenger supermultiplet couples at tree level to the spurion and through 
the MSSM gauge fields (in this diagram the B boson) with the MSSM. When the messengers are 
integrated out from the action, this diagram contributes to the soft mass of the up squark. 

Minimal anomaly mediation. The advantage of AMSB over SUGRA is that it natu- 
rally conserves flavour. However, pure anomaly mediation leads to tachyonic sleptons, i.e. 
their squared soft masses become negative. This would cause them to acquire non-zero 
VEVs and break the electromagnetic gauge symmetry. The minimal AMSB (mAMSB) 
model uses a phenomenological approach to tackle this problem: a universal, non-anomaly- 
mediated contribution TOq is added to the soft squared scalar masses (3.3b) at the scale 
M AMSB . The origin of these terms may be for example additional fields in the bulk, but in 
the mAMSB model mo is simply considered as a parameter of the model. 

3.7 General gauge mediation 

Several models of Gauge-Mediated Supersymmetry Breaking (GMSB) have been proposed 
in the literature (see [32] for a review). Many of these models include a field X, called 
the spurion, that acquires a supersymmetry-breaking VEV, and a set of weakly coupled 
fields that are charged under the MSSM. The latter are called messenger fields since they 
communicate supersymmetry breaking to the MSSM fields: they interact at tree level with 
the spurion and through the MSSM gauge fields with the MSSM (see figure 3). 

Recently, the framework of General Gauge Mediation (GGM) [33, 34] has been pro- 
posed to unify all earlier descriptions of GMSB. It describes the effects of an arbitrary 
hidden sector on the MSSM. It starts from the following definition of gauge mediation: in 
the limit of vanishing gauge couplings, the theory decouples into the MSSM and a separate, 
supersymmetry-breaking hidden sector. For example, the setup described above fits into 
this definition by taking the messenger and spurion fields as the hidden sector. 

In the GGM framework, all MSSM soft terms can be described in terms of a small 
number of correlation functions involving hidden-sector currents. Essentially, the GGM 
framework parametrises the effects of the hidden sector on the MSSM. By constructing the 
effective Lagrangian, the following soft-mass formulae are found: 4 



M a = g 2 a B, 



'a 



(a = 1,2, 



3) 



(3.6a) 



3 




(3.6b) 



a=l 



4 For future convenience, a factor M (the messenger scale) has been absorbed into the definition of 
the B a , cf. [9]. 
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with 

Ci(!) = |sf, (3.7a) 

aw-/* for *< = Q. £ : ff ..^. (3.7b, 

[0 for $j = u R , d R , e R , 

[I for $j = Q, mb, dp, 

C 3 (i) = < 3 ~ _ ( 3 - 7c ) 

[0 for $i = L,e Rl H u ,H d . 

Here £? a , £ and A a are expressions involving the hidden-sector current correlation functions; 
Yi is the hypercharge of the scalar field and C a (i) is the quadratic Casimir of the 
representation of <&j under the gauge group labeled by a. Usually a Z2 symmetry of the 
hidden sector is assumed in order to forbid the term containing £, since it would lead to 
tachyonic sleptons. The above conditions are the matching conditions at the messenger 
scale M where we integrate out the hidden sector. The seven numbers (,A a ,B a contain 
information on the hidden sector, but are regarded as parameters of the low-energy EFT 
that we call the MSSM. 

The GGM framework does not allow for additional interactions that could generate 
\jl and b radiatively; that would require interactions between the MSSM and the hidden 
sector that remain in the limit of vanishing gauge couplings. The framework would have to 
be extended to allow for such couplings. To parametrise the effects of such an extension, 
additional contributions 5 u ,8d to m 2 Hu ,m 2 Hd are often added. 

Minimal gauge mediation. Minimal gauge mediation (MGM) is a GGM model that 
is restricted to a subset of the GGM parameter space, defined by the constraints A a = A, 
B a = B and A = IB 2 . The term corresponding to £ is taken to be zero. Additional 
non-gauge contributions 5 U , 84 are added to the soft Higgs masses. Then the expressions 
for the soft masses become 

M a = g 2 a B a (a = 1,2, 3), (3.8a) 



m 2 i =2B 2 Y J 9tC a {i), (3.8b) 

a=l 

3 

m 2 Hu = 2B 2 dtCa(Hu) + S u , (3.8c) 

a=l 

3 

m 2 Hd = 2B 2 Y J 9tC a {H d ) + 6 d , (3.8d) 



a=l 



where this time mf denote only the squared masses of the squarks and sleptons. 
3.8 Mirage mediation 

Rather than restricting oneself to one of the three known mechanisms for radiative super- 
symmetry breaking (gravity, anomaly or gauge mediation), one could solve the problems of 
particular models by choosing two (or more) mechanisms and combining the best of both 
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worlds. For example, one might tackle the tachyonic slepton problem of anomaly mediation 
by combining it with gauge mediation (see e.g. [35]). 

Mirage mediation [36] is one such scenario in which gravity-mediated and anomaly- 
mediated soft terms have comparable contributions. In this scenario, the gravity-mediated 
terms are suppressed by a relative factor log (Mpi/m 3/ / 2 ), which is numerically of the order 
of a loop factor. This results in mirage unification: the gaugino and scalar masses unify 
at a scale far below the scale where the soft masses are generated. This mirage messenger 
scale does not correspond to any physical threshold, hence the name. 

This class of phenomenological models are based on a class of string models with 
stabilised moduli, called the KKLT construction. It solves the tachyonic slepton problem 
that arises in pure anomaly mediation and has reduced low-energy fine-tuning [37]. 

3.9 One-loop RGIs for the MSSM 

It should be noted that if we have a set of RGIs, then any function of those RGIs will 
also be RG invariant. Therefore, in order to find all RGIs, one should look for a maximal 
set of independent RGIs, i.e. invariants that cannot be expressed in terms of each other. 5 
Recently a complete list of independent one-loop RGIs for the MSSM was derived in [7, 8]: 

D Bl , = 2 (m^ - m|J - m| + m\ - m| + m|, (3.9a) 
D Ll3 = 2 (mi - m|J - m| + m|, (3.9b) 
D X1 = 3 (3m| - 2 (ml - m|J - m\^ - m|, (3.9c) 

2_ 

Qi 
10 



Dv, , cr = rn ~ — 2ml + ml — m% + ml 



m~ 



... , 2ml + ml — m% + ml + m?r — m?r ), (3.9d) 

13 V Qa «s d 3 l 3 ' e 3 tin H d j ' v ; 

D Z = 3 (m| - m|) + 2 (m| 3 - m 2 H J , (3.9e) 



^ = 72 ( m L " m H d + H 
"1 \ gen 



m~ - 2m| + ml - m| + m|J J , (3.9f) 

^ gen / 

lB a = ^ (a = 1,2,3), (3.9g) 

9 a 

33/2 2 2 



/ Ml =Mt- T [mA^ - - mfj , (3.9h) 

= Mi + 1 (9 (ml - ml) + IGm^ - m|) , (3.9i) 

/m 3 = M| - ^ ( 5m | + ^ _ m |) , ( 3 .9j) 
33 

/ g2 =5r 2 -y<?2" 2 , (3.9k) 

^3=9r 2 + y5 3 " 2 5 (3-91) 



3 It is tempting to call this a 'basis of RGIs', as in [8]. Note however that it is not the same as a basis 
of a vector space. One should keep in mind that once we have found such a set, we are not restricted to 
making linear combinations of them, but can also take products, quotients and so on. 
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9 256/3 \ V 61 
929 3 \ 



(3.9m) 



27 „ 21 10 , 256 _ 9 , r 73 , r 

i 4 = ^ f ^ b A T M 3 M 2 H Mi. (3.9n) 

/i 61 61 61 183 61 2013 v ; 

The sum in (3.9f) runs over the three sfermion generations. The notation has been taken 
over from [8] (equations (3.9a)-(3.91)) and [7] (equations (3.9m)-(3.9n)). These one-loop 
RGIs have been derived under the pMSSM assumptions, using the relevant /^-functions 
listed in appendix A. An alternative, systematic way of deriving this complete set can be 
found in appendix B. 

Any other (one-loop) RGI we can think of can be written in terms of those in (3.9). 
For example, the RGI in the example from section 2.4.3 can be written as 

. _9 48 n 28 a 

I g = 4 5l - —9 2 + y 9 3 

= n J » + ii /s3 - (3 - 10) 

The two RGIs in (3.9m)-(3.9n) will not be relevant to our analysis. We will explain why 
in the next subsection. 

3.10 RGIs in the literature 

As we have seen, RGIs provide a new tool to test predictions about high-scale physics, such 
as gauge-coupling unification. The trick is to find sum rules for high-scale physics that can 
be written in terms of RGIs. In the literature, several such sum rules can be found. 

Consider for example minimal gauge mediation, which has been studied in the context 
of RGIs in [9]. If one inserts the spectrum (3.8) at the messenger scale into the RGI 
expressions from (3.9a)-(3.91), one immediately finds that Db 13 = D^ 13 = D Xl = 0. In 
terms of the model parameters of MGM, the non-vanishing RGIs are 

D Y 13H = -^(S u -5 d ), (3.11a) 

D z = -28 d , (3.11b) 

I Ya = gf{M) (6 U - 8 d ) , (3.11c) 

I Ba =B (a = 1,2,3), (3.11d) 
38 



OQ 

I Ml = -g\{M)B\ (3.11e) 



I M2 = 2gl(M)B 2 , (3.11f) 

Im 3 = -2gl(M)B 2 , (3.11g) 
33 

I g2 = g -\M)-- g2 -\M\ (3.11h) 

I g3 = g -\M) + ^g^(M), (3-lli) 

where M is the messenger scale. This gives us eleven equations in terms of six unknowns 
(5 U , 5d, B, <7i(M), g2 (M), g^ (M)). We can trade each unknown for an equation, i.e. for 
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each parameter we use one of the above equations to express it in terms of RGIs only. 
Since we have more independent equations than unknowns, we can substitute the resulting 
six expressions into the remaining five equations to obtain five sum rules in terms of RGIs 
only. For example, using equation (3. lid) we can eliminate the model parameter B from 
the remaining equations. Then equations (3.11e)-(3.11g) can be used to eliminate the gauge 
couplings at the messenger scale. We can get the value of Sd from equation (3.11b), and 
then (3.11a) gives the value of S u . After substituting the resulting six expressions into the 
five remaining equations, we are left with the following sum rules: 




13 / 38 

= I Bl -I B2 , (3. lid) 

(3.11d) 

I g2 , (3.11h) 

I g3 , (3-lli) 

To summarise, we have chosen a specific supersymmetry-breaking model and expressed 
the RGIs in terms of model parameters. Since we ended up with more equations than 
unknowns, we could construct eight sum rules in terms of RGIs only: three from vanishing 
RGIs and five by eliminating the model parameters. If any of these sum rules are violated, 
MGM is not consistent with experimental data. This test can be performed at any energy 
scale where the MSSM is valid, which implies that we don't have to know the messenger 
scale M in order to rule out the MGM model. 

Now we can see why the two RGIs in (3.9m)-(3.9n) are not useful. Suppose we wish 
to test a specific breaking model. Let us denote the values of B = b/fi and \i at the 
new-physics threshold in this model as B t hr and /ithr respectively. Now we apply the above 
procedure to this model: we express _B t hr and /xthr m terms of RGIs and the other couplings 
at the high scale; then we can insert these expressions into the remaining equations. But 
since B and [i both appear in only one independent RGI, there are no equations to insert 
these expressions into! In the above example, it was possible to combine all RGIs into 
sum rules because each running parameter appeared in more than one RGI. Since B and 
\x do not, their corresponding RGIs become useless to our analysis. Hence we will have to 
restrict ourselves to the RGIs that do not contain these running parameters. As will be 
shown in appendix B, it is not possible to construct RGIs that contain the Yukawa and soft 
trilinear couplings without using B and /u. That is why we will only use RGIs constructed 
out of soft masses and/or gauge couplings, i.e. the RGIs listed in (3.9a)-(3.91). 

3.10.1 Using RGIs effectively 

In studies of RGIs such as [8-10], a certain breaking mechanism is usually presupposed. 
Then one constructs sum rules that are tailor-made for that breaking mechanism. For 
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example, the sum rules constructed above all provide a test for consistency of MGM with 
experimental data. However, some of these sum rules will also hold for other breaking 
mechanisms. It is not always clear to what extent the validity of the sum rules depends 
on the unique features of the breaking mechanism under study. For example, in minimal 
gauge mediation the quantity M a jg\ unifies at the messenger scale; this follows directly 
from the matching condition (3.8a). However, in mSUGRA this quantity also unifies, but 
for a different reason: it is the consequence of the assumption of gauge-coupling unification 
and gaugino-mass unification at the same energy scale! Hence, the sum rules that test 
this unification property cannot be used to confirm that either of these specific models 
corresponds to reality. They can only provide consistency checks that should be satisfied 
if any of these models are realised in Nature. 

Therefore, we will look for RG-invariant sum rules using a different, more model- 
independent approach. We will not presume any spectrum specific to a certain breaking 
mechanism. Instead, we will search for sum rules that test properties that are common in 
supersymmetry-breaking models (e.g. M a jg\ unification). Then any breaking model, be it 
an existing one such as those described in sections 3.5-3.8 or a new one contrived in the 
future, can be tested directly if it predicts any of these properties. For example, if the 
sum rules for M a jg\ unification are not satisfied by experimental data, then models that 
predict this property (mSUGRA and MGM, but not necessarily GGM) are falsified. But 
also anyone who would concoct a new model that has this property, would have to go back 
to the drawing board at once. In the next section, we will look for common properties to 
test and find sum rules for them. 

4 Results 

Supersymmetry-breaking models predict relations between the running parameters as a 
result of matching conditions at the new-physics threshold. These relations mostly involve 
the unification of certain parameters. Therefore we will construct sum rules for the following 
scenarios. 

Scenario 1: gauge-coupling unification. As can be seen from figure 2, the MSSM 
may be consistent with gauge-coupling unification, depending on the values of the sparticle 
thresholds. The hypothesis that the gauge couplings unify is often made in supersymmetry- 
breaking models, for example in mSUGRA. Therefore it will be important to determine 
whether gauge-coupling unification occurs in Nature. We will call the special case where 
9i= 92 = 93 = 9 scenario 1. 

Scenario 2: gaugino-mass unification. In mSUGRA, the gaugino masses are assumed 
to unify. Since this model is widely used, it is useful to check whether gaugino-mass unifi- 
cation occurs in Nature. We will call the case where M\ = M2 = M3 = M\/2 scenario 2. 

Scenario 3: Unification of M a /g^. As we mentioned in section 3, the quantities 
Mi/gf, Mijg\ and M^jg\ may unify for different reasons. It could be the consequence of 
gaugino-mass and gauge-coupling unification at the same scale (as in mSUGRA) or it may 
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be the result of the gaugino-mass matching conditions (as in MGM). Therefore we will also 
test this property. We will call the case where M\jg\ = Mijg\ = -^3/53 = C scenario 3. 



Scenario 4: flavour-universality of high-scale sfermion masses. In many theories, 
the sfermion masses are assumed to be flavour-universal, i.e. the first (and therefore also 



This hypothesis is motivated by the need to suppress FCNC amplitudes. Flavour-universality 
may be a consequence of flavour symmetries (as postulated in mSUGRA) or of the flavour- 
blindness of the interactions that mediate supersymmetry breaking (as in GGM). Since 
this property occurs in many models, we will test flavour-universality of the soft masses. 

Scenario 5: scalar-mass unification. Unification of the soft scalar masses, which 
occurs in mSUGRA, is very predictive: many matching conditions depend on a single pa- 
rameter mo, which allows us to construct multiple sum rules. Sometimes non- universality 
of the soft Higgs masses is assumed, because suppression of FCNC amplitudes does not 
require them to be universal with the squark and slepton masses. Therefore we will dis- 
tinguish between two cases. We will refer to the case where w? H ^ and mr H have additional 
non-universal contributions 5 U and 5^ respectively as scenario 5a. The special case of 
universal scalar masses, i.e. 5 U = 5d = 0, will be denoted as scenario 5b. 

Mixed scenarios: multiple unifications at one scale. It is possible that several 
of these unification properties will turn out to be consistent with experimental data. In 
that case, we could test whether these unifications occur at the same energy scale. At 
first sight, it may seem strange to consider the possibility of two kinds of unifications 
at different energy scales. In supersymmetry-breaking models such unifications usually 
occur at a threshold where new physics enters the theory. Thus even if the MSSM were 
consistent with two kinds of unifications, the RG trajectories of the running parameters 
could be deflected from the MSSM trajectories after the first threshold, spoiling the second 
unification. However, recall that in mirage mediation (see section 3.8) the scale where the 
soft masses unify is lower than the scale at which the soft masses are generated. Thus it is 
possible that the unification scale does not correspond to any physical threshold. Therefore 
we will separately check whether multiple unifications occur at the same scale. 

4.1 Sum rules 

In this section we construct sum rules that test the scenarios described above. 

Scenario 1. From equations (3.9a)~(3.91) and table 3 we can see that we have 14 equa- 
tions with 16 unknowns, so at first sight we expect to find no sum rules. However, equations 
(3.9k) and (3.91) form a subset of two equations with only one unknown g. Hence we can 
make one sum rule: 




^92 + 4-^33 — 0- 
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Scenario 2. In this scenario we also have 14 equations with 16 unknowns. However, 
equations (3.9g), (3.9k) and (3.91) form a subset of five equations with four unknowns. 
This allows us to construct one sum rule: 

33 \ f 11 \ 

I Bl ~ y 1 ^) h, = ( + -jIbA I 92 - (4.2) 

Scenario 3. Again we have 14 equations with 16 unknowns. This time equations (3.9g) 
form a subset of three equations with one unknown. This yields two sum rules: 

I Bi =Ib 2 , (4.3) 
I Bi =Ib 3 . (4.4) 

Scenario 4. We can see directly from table 3 that this scenario yields two sum rules in 
the form of vanishing RGIs: 

D Bl3 = 0, (4.5) 
D Ll3 = 0. (4.6) 

The remaining 12 equations contain 13 unknowns and no subset of them contains less 
unknowns than equations. Hence, we cannot construct any other sum rules. 

Scenario 5a. In this scenario we also have two vanishing RGIs, because unified scalar 
masses imply flavour-universal scalar masses. The remaining RGIs yield twelve equations 
with nine unknowns. Hence, we can construct three sum rules that are specific to this 
scenario: 

/ 33 \-l/2 33 / i \-i/2 

I g2 = I Bl \ I Mi - ^ D xi) - Y Ib2 \ Im2 ~ 8 Dxi J ' ( ' 4 ' 7 ' ) 
/ 33 \-V2 11 / 3 \-V2 

^=^(^-^1) +-5-^(^3 + ^x1) > ( 4 - 8 ) 

I 33" 13 

= I Ya ^h h ~ -^D X1 + -I Bl D Yl3H . (4.9) 

Furthermore, non-universality of the Higgs masses can be tested directly because we can 
extract 5 U and from the RGIs: 

5 d = -\d z + 0, (4.10) 



13 

5u = ~W DY ^ H ~ 2 Dz * °' 
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Invariant 


Scenario 1 
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-^92 
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^92 


^92 
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^93 
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Table 3. Values of the MSSM RGIs in the unification scenarios 1-5. For each scenario, RGIs that contribute to sum rules that are specific to that 
scenario arc listed in boldface. If an RGI is simplified with respect to its definition but is not used for a sum rule, its name with the appropriate 
substitutions is listed. For the substitution (m~ — > m~) it is implied that / = Q, u, d, L, e and i = 1, 2, 3. If an RGI does not simplify at all, only 
its name is listed. 



Scenario 5b. We can see directly from table 3 that this scenario yields three new sum 
rules in the form of vanishing RGIs: 



D Yrm = 0, (4.12) 
D z = 0, (4.13) 
I Ya = 0. (4.14) 

The remaining nine equations contain seven unknowns, which allows us to construct two 
more sum rules: 

/ 33 \~ 1/2 33 / 1 \ _1 /2 

^=^k-^a) ~Y lBa \ lM2 ~8 Dxx ) ' (415) 

/ 33 \ - 1 / 2 11 / 3 \ ~ 1/2 

I 93 = I Bl (l Ml - -D Xl j + -I Ba [I M3 + -D Xl j . (4.16) 

Note that these two sum rules also hold for scenario 5a. This makes sense, because sce- 
nario 5b is a special case of scenario 5a, so the sum rules (4.7)-(4.9) will also hold for 
scenario 5b. However, sum rule (4.9) has become redundant because it is automatically 
satisfied if (4.12) and (4.14) hold. 

Mixed scenario 123. Note that if any two of the scenarios 1, 2 and 3 hold, the hypothesis 
that the corresponding unifications occur at the same scale is equivalent to the hypothesis 
that the third scenario also holds. For example, if the gauge couplings and the gaugino 
masses unify, these unifications occur at the same scale if and only if M a /g 2 unification 
occurs. Hence, for simultaneous unification at least the sum rules (4.1)-(4.4) should hold. 
Note that (4.2) is automatically satisfied if the other three sum rules hold, so we have three 
independent sum rules for this mixed scenario. 

From table 3 we conclude that if these three unifications occur simultaneously we have 
14 equations with 14 unknowns, which gives no sum rules at first sight. However, equations 
(3.9g), (3.9k) and (3.91) constitute five equations with two unknowns, giving us three sum 
rules. We have already found three, so there are no new sum rules. 

Mixed scenario 15. If both scenario 1 and scenario 5a hold, at least the sum rules (4.1) 
and (4.5)-(4.9) hold, as well as the inequalities (4.10)-(4.11). From table 3 we conclude 
that if these unifications occur simultaneously we have two sum rules from vanishing RGIs 
as well as twelve equations with seven unknowns. Thus we can construct five additional 
sum rules. We already found six of them, so there is one new sum rule: 



13 

lY a = ^I 92 D Yl . AH . (4.17) 



If both scenario 1 and scenario 5b hold, we have Dz = Iy a = Dy-izh = 0- ^ n that case this 
new sum rule becomes redundant. 

Mixed scenario 25. If both scenario 2 and scenario 5a hold, at least the sum rules (4.2) 
and (4.5)-(4.9) hold, as well as the inequalities (4.10)-(4.11). From table 3 we conclude 
that if these unifications occur simultaneously we have two sum rules from vanishing RGIs 
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as well as twelve equations with seven unknowns. Hence we can construct five additional 
sum rules. We already found six of them, so there is one new sum rule: 

81 56 

hh ~ ^hi 2 + 25/M3 = (4.18) 
Note that we get the same additional sum rule if we choose scenario 5b instead of 5a. 

Mixed scenario 125. If all unifications occur simultaneously, then we immediately know 
that the following sum rules should hold (we assume scenario 5a for the moment): 

• Sum rule (4.1) for scenario 1. 

• Sum rules (4.3)-(4.4) for scenario 3. 

• Sum rule (4.2) for scenario 2. However, as we have mentioned before, this one is 
automatically satisfied if the sum rules for scenarios 1 and 3 are satisfied. Hence this 
one is redundant. 

• Sum rules (4.5)-(4.6) for scenario 4, which is implied by scenario 5. 

• Sum rules (4.7)-(4.9) for scenario 5, as well as inequalities (4.10)-(4.11). 

• Sum rule (4.17) for the simultaneity of scenarios 1 and 5. 

• Sum rule (4.18) for the simultaneity of scenarios 2 and 5. However, this one has 
become redundant: it can be retrieved by combining equations (4.3)-(4.4), (4.7)-(4.9) 
and (4.17). We could have expected this: if we have consistency with scenarios 1, 2 
and 5, and we have established simultaneity of both scenarios 1 and 2 and scenarios 1 
and 5, then it follows automatically that we have simultaneity of scenarios 2 and 5. 

This adds up to nine independent sum rules. From table 3 we conclude that if all unifica- 
tions occur simultaneously we have two sum rules from vanishing RGIs as well as twelve 
equations with five unknowns. Hence, we should get nine sum rules, so there are no new 
ones. Note that if we had taken scenario 5b instead of 5a, we would not get any additional 
new sum rules either. In that case the usual sum rules Dz = Iy a = Dy 13H = would hold 
and equation (4.17) would become redundant. 

4.1.1 Sum rules summary 

All scenarios discussed above and their corresponding sum rules have been summarised 
in figure 4. Related scenarios have been connected: if one starts at a given scenario, one 
should follow the arrows downwards to arrive at the underlying hypotheses. When we have 
determined the values of the RGIs from experimental data, we can test whether the listed 
scenarios are consistent with the data. One should proceed as follows: to test a hypothesis, 
check the validity of the sum rules in the corresponding box. Then check the validity of 
the sum rules in all boxes one encounters by following the arrows all the way down. If all 
these sum rules are satisfied, the hypothesis is consistent with the experimental data (as 
far as our sum rules are concerned). 
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Figure 4. Scheme for testing hypotheses about the spectrum at the new-physics threshold. For a 
given scenario, the arrows point towards its underlying hypotheses. To test a specific hypothesis, 
check whether the corresponding sum rules arc satisfied. Then follow the arrows downwards all the 
way to the bottom and for each sum rule along the way, check whether it is satisfied. 
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4.2 Model-specific sum rules 

Until now we have only considered hypotheses concerned with relations between the run- 
ning parameters of the (p)MSSM. These hypotheses do not refer to any model-specific 
parameters. However, we can find additional sum rules for certain models because the 
soft masses are related by only a few parameters. For example, in MGM the gaugino and 
sfermion masses are determined by the gauge couplings and a single parameter B; see 
equations (3.8a)-(3.8b). Furthermore, the question whether the messenger scale equals the 
gauge-coupling unification scale only makes sense if we consider gauge-mediation models. 
Therefore we consider model-specific sum rules separately in this section. 

In the following, we will look for sum rules for GGM and AMSB that do not follow 
from the general hypotheses we have discussed above. Such sum rules will generically be 
referred to as "model-specific sum rules". The sum rules for GGM can also be found in [9]. 
To our knowledge no sum rules for AMSB models have been presented in the literature. We 
will not discuss mSUGRA, because for our purposes the mSUGRA spectrum is completely 
characterised by simultaneous scalar-mass, gaugino-mass and gauge-coupling unification. 
So, there are no mSUGRA-specific sum rules. 

4.2.1 General gauge mediation 

In section 3.7 the RG boundary conditions for GGM are given in terms of 11 model param- 
eters at the messenger scale M: 8 U ^ and B a , A a , g a (M) for a = 1, 2, 3 (we take £ = 0). If we 
insert these RG boundary conditions into the RGIs, we immediately find Db 13 = Dl 13 = 0, 
as expected from flavour-universality. In addition we get one model-specific sum rule: 

D X1 = 0. (4.19) 

The remaining RGIs have the following values: 

DY 13H = -^(S u -6 d ), (4.20a) 

D z = -25 d , (4.20b) 

Ir a = 9i 2 (S u - 5 d ) , (4.20c) 

I Ba =B a (a = 1,2,3), (4.20d) 

Im* = A [b\ + , (4.20e) 
1 

2' 

I M3 = gj (b* - , (4.20g) 
33 

I 9 2=9i 2 --^92 2 , (4-20h) 

= 9i 2 + Y 9 * 2 ' (4 - 20i) 

where the gauge couplings are understood to be evaluated at the messenger scale. This 
amounts to eleven equations with eleven unknowns, hence no additional sum rules can 



Im 2 = gt ( B\ + \A 2 ) , (4.20f) 
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be constructed. Note that we can again verify non-universality in the Higgs sector using 
equations (4.10)-(4.11). 



lY a = ^I 92 D Yl3H . (4.21) 



Gauge-coupling unification at the messenger scale. If both gauge-coupling unifi- 
cation and GGM are compatible with experimental data, we may ask ourselves if the mes- 
senger scale equals the scale of gauge-coupling unification. If we insert g\ = g 2 = 93 = 9 
into equations (4.20), we get eleven equations with nine unknowns. Hence we can make 
two more sum rules. This includes equation (4.1) for gauge-coupling unification. Hence, 
there is only one additional model-specific sum rule: 

13 

56 J 

Note that this sum rule happens to be identical to equation (4.17), which was constructed 
for testing a completely different concept: simultaneous scalar-mass and gauge-coupling 
unification (i.e. scenario 15). Therefore we have marked this sum rule as being model- 
specific. Also note that this sum rule becomes redundant in the case of universal Higgs 
masses, since Iy a = Dy 13H = in that case. 

Minimal gauge mediation. Recall that MGM is a GGM model restricted to a subset 
of the GGM parameter space defined by A a = A, B a = B and A = 2B 2 . Inserting this 
into the RGI values (4.20) of GGM, we find eleven non-vanishing RGIs that depend on six 
parameters. Hence we can construct five sum rules. These include (4.3)-(4.4) for M a /g 2 
unification. There are three additional model-specific sum rules:® 



13 

+ > ( 4 - 22 ) 



38 33^ / 2 
5/mi 5 1 V Im 2 




° = Wt7 —tBlJ-f T 92, (4-23) 



<424) 

Note that (4.22) becomes redundant in the case of universal Higgs masses. 
4.2.2 Anomaly mediation. 

The RG boundary conditions for AMSB are given in section 3.6 in terms of 4 model 
parameters at the scale M AMSB where supersymmetry breaking occurs: m 3 / 2 an d 9a(M AMSB ) 
for a = 1,2,3. If we insert these RG boundary conditions into the RGIs, we immediately 
find nine model-specific sum rules: 

D Bl3 = D Ll3 = D X1 = D Yl3H = D z = ly a = hh = hi 2 = hi 3 = 0. (4.25) 



6 Note that in MGM, we can safely divide by lM a '- if one of the Jjvj a vanished, then B — and the 
gaugino masses would vanish at the messenger scale. Their /3-functions, being proportional to the gaugino 
masses, would vanish as well. Then at one-loop order gauginos would be massless at all scales (only through 
two- loop effects the masses will be non- vanishing) . In that case we would have observed them already. Thus 
the Im„ cannot vanish. 
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Note that Db 13 and Dl 13 vanish although the sfermion masses are not flavour-universal! 
The non-vanishing RGIs have the values: 

T 33 m3 / 2 (A Oft \ 

I Bl = ~— 2 , (4.26a) 
Ib 2 = (4.26b) 

Ib ^ = -^B^ (4 - 26c) 

^=5r 2 -y5 2 - 2 , (4.26d) 
^3=5r 2 + y^ 2 , (4-26e) 

where the gauge couplings should be evaluated at the scale of supersymmetry breaking. 
This amounts to five equations with four unknowns, but we can do better: equations 
(4.26a)-(4.26c) constitute three equations with one unknown. This yields another two 
model- specific sum rules: 

33 

= I Bi - I B2 , (4.27) 
5 

= /B 1 +y/B 3 - (4-28) 

Gauge-coupling unification at the scale of supersymmetry breaking. If both 
AMSB and gauge-coupling unification turn out to be consistent with experimental data, we 
may ask ourselves whether supersymmetry breaking occurs at the scale of gauge-coupling 
unification. In that case we should insert g\ = 92 = 92 = 9 into (4.26). But this will 
only affect equations (4.26d) and (4.26e), which we have not used to make the above sum 
rules. This amounts to two equations with only one parameter, so we get one more sum 
rule. This must be the sum rule (4.1) for gauge-coupling unification, hence there are no 
sum rules that specifically test whether gauge-coupling unification occurs at the scale of 
supersymmetry breaking. 

Minimal anomaly mediation. Recall that in minimal AMSB, a universal additional 
term m\ is added to the soft scalar masses. If we insert this into the RGI expressions, we 
immediately find five model-specific sum rules: 

D Bi:i = D Ll3 = D Yl3H =D Z = I Ya = 0. (4.29) 
The non-vanishing RGIs have the values: 

D X1 = 5ml, (4.30a) 

I Bl = ^T7T2, (4.30b) 

5 lbTT z 

I B2 = j^, (4.30c) 

!b* = (4.30d) 
33 

Imi = -s-mg, (4.30e) 
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I Ma = ^ml (4.30f) 

J Ms = -^rnl (4.30g) 

I g2 =gf--g^, (4.30h) 

^ 3 =5r 2 + y5 3 " 2 , (4.3(H) 

where again the gauge couplings should be evaluated at the scale of supersymmetry break- 
ing. This adds up to nine equations with five unknowns, so we expect to find four additional 
sum rules. However, if we leave out equations (4.30h)-(4.30i), we are left with seven equa- 
tions with only two unknowns. This yields another five model- specific sum rules: 



= I Bl -—I B2 , (4.31) 
5 

= ^ + ^3, (4.32) 
40 

= D X1 - —I Ml , (4.33) 

= D Xl -8I Ah , (4.34) 
1 6 

= D xl + -I Ms . (4.35) 



Here equations (4.31)-(4.32) also hold for AMSB. Equations (4.33)-(4.35) are automatically 
satisfied in AMSB because D Xl and Iu a vanish. 

4.3 Discussion 

In this section, we have found (a) sum rules that test general properties of the RG boundary 
conditions and (b) sum rules that test the consistency of specific model spectra. Comparing 
both sets of sum rules will help us to determine how good the sum rules are at distinguishing 
between several properties and model spectra. 

Sum-rule ambiguities and how to eliminate them. In the sum rules we observe the 
following ambiguities: 

• If the sum rules (4.1) for gauge-coupling unification and (4.3)-(4.4) for M a jg\ uni- 
fication are both satisfied, then the sum rule (4.2) for gaugino-mass unification is 
automatically satisfied. But gaugino-mass unification is implied by gauge-coupling 
unification and M a jg\ unification only if both unifications occur at the same scale! 
Hence, if (4.1), (4.3) and (4.4) are satisfied by experimental data, then we cannot 
determine unambiguously whether the gaugino masses unify. At this point, we should 
use the bottom-up method to examine the running of the parameters. Then we could 
see whether the unification scales are the same. 

• Equation (4.17) checks whether scalar masses and gauge couplings unify at the same 
scale. Equation (4.21) checks whether the gauge couplings unify at the messenger 
scale in GGM. Yet these sum rules happen to be the same. However, this does not 
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mean we cannot distinguish between these two scenarios. The former scenario also 
requires that the sum rules (4.7)-(4.9) for scalar-mass unification are valid. In the 
latter scenario, these sum rules are not satisfied. Thus the double role of (4.17) poses 
no problem. 

• In AMSB and mAMSB, the sum rules (4.5)-(4.6) for flavour-universality are satisfied, 
although the sfermion masses in these models are clearly non-universal. Fortunately, 
(m)AMSB has a lot more sum rules, which could help discern these models from 
flavour-universal ones. For example, the vanishing of Dy 13H , Dz and Iy a is typical 
for (m)AMSB. Equations (4.7)— (4.8) then help us discern (m)AMSB from scalar-mass 
unification with universal Higgs masses. Again, satisfying a single sum rule may be 
ambiguous, but other sum rules eliminate this ambiguity. 

• Because Dy 13H and Iy a vanish in (m)AMSB, the sum rule (4.17) for simultaneous 
gauge-coupling and scalar-mass unification is automatically satisfied. However, the 
sum rules for gauge-coupling unification and scalar-mass unification again help us 
distinguish between both scenarios. 

• The vanishing of Iu a in AMSB and the sum rules (4.33)~(4.35) of mAMSB both imply 
that the sum rule (4.18) for simultaneous scalar- mass and gaugino-mass unification 
is satisfied. However, the sum rules for scalar-mass unification and gaugino-mass 
unification help us distinguish between both scenarios. 

If we only consider the spectrum properties and breaking mechanisms that we dis- 
cussed in this section, our sum rules work surprisingly well. Many of the sum rules are not 
unambiguous by themselves, but in most cases the other sum rules remove the ambiguity. 
Only when the data are consistent with both gauge-coupling unification and M a jg\ unifi- 
cation, we have to resort to other methods (such as the bottom-up method) to determine 
whether the gaugino masses also unify (or equivalently, whether both unifications occur at 
the same scale). 

Of course, it is possible that a new supersymmetry-breaking model is concocted in the 
future, and that some of its corresponding sum rules introduce similar ambiguities. These 
may or may not be resolved by other sum rules. Therefore, we should keep in mind that 
if the sum rules of a model or hypothesis are satisfied, this is not a confirmation that this 
model or hypothesis is correct. The true power of our sum rules is their falsifying power: 
the failure to satisfy just one sum rule implies that the corresponding hypothesis or model is 
incorrect. This falsifying power depends on how precisely the parameters entering the RGIs 
have been measured: the smaller the errors are, the more powerful the rejection becomes. 

Now that we have an idea of the quality of the RGI sum rules, we can finally examine the 
advantages and limitations of the RGI method. 

Advantages of RGIs 

• The RGI method requires less input than the other methods we have discussed. We 
only need the values of all soft masses and gauge couplings at one scale. These are 
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sufficient to reconstruct the values of the RGIs in equations (3.9a)-(3.91). In contrast 
to the bottom- up method, we do not need the values of the Yukawa couplings and 
soft trilinear couplings because we could not use them anyway. Also, the value of the 
new-physics threshold does not have to be known. 

• The RGI method is very simple: it is entirely algebraical and does not require the 
numerical integration of renormalisation group equations. Therefore it avoids the 
complicated propagation of errors between the collider scale and the new-physics 
threshold. Also, it is not as time-consuming as the top-down method. 

• If no unification of any sort is observed, the bottom-up method provides us with no 
clues about new physics, whereas the RGI method may still be helpful. Suppose 
for example that non-minimal General Gauge Mediation is a good model of Nature, 
so that no unification occurs. Then the bottom- up method does not provide any 
hints on where to look for new-physics thresholds. We could check whether the 
sfermion masses match the pattern (3.6b) at any scale, but that would take a lot of 
time. However, if we apply the RGI method we only need to check the three sum 
rules (4.5), (4.6) and (4.19) to see whether this bottom-up scan is worth the effort. 

• As long as just a few of the relevant soft masses have been measured experimentally, 
the sum rules can be exploited as a fast means of identifying theoretically interest- 
ing regions in the remaining parameter space (e.g. regions with specific unifications). 
Here we assume that the particles corresponding to the masses have been identified 
correctly. The precise procedure for reconstructing the MSSM parameters from ex- 
perimental data is beyond the scope of our study and will be investigated elsewhere. 

Limitations and challenges of the RGI method 

• As we mentioned before, RG invariance only holds up to a certain loop level. The 
RGIs in (3.9) have been determined using the one-loop RG equations. Higher-order 
loop effects will certainly spoil RG invariance. We could of course try to find RGIs 
for the MSSM at a higher loop order. But already at the two-loop level the RG 
equations for the MSSM (see e.g. [38]) are too complicated to retain the simplicity 
of this method, if it is possible to find RGIs at all. 

However, the relevant question is to what extent we should worry about this approx- 
imate RG invariance. It has been demonstrated in [8] for the MSSM that even if we 
have determined all soft masses at the collider scale with a universal experimental 
uncertainty as low as 1%, two- loop contributions to the RGIs are smaller than the 
experimental errors in the RGIs. Thus for all practical purposes we can safely treat 
the one-loop RGIs as true invariants. It should be noted that if we were to apply 
the RGI method to an EFT other than the MSSM, we should again study the size of 
2-loop contributions to the RGIs before we can neglect them. 

• It may seem like the RGI method magically reduces the uncertainties of the running 
parameters, compared to RG-evolved parameters. However, we have paid a price for 
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this reduction, namely information. We can directly see this from equations (3.9a)- 
(3.91): we started with 18 running parameters (12 scalar masses, 3 gaugino masses 
and 3 gauge couplings) and have reduced them to only 14 invariants. 

We can easily understand why we have to give up information to gain smaller errors. 
Consider for example the RG equations for m?~ and m 2 ~ (see appendix A for the 
definitions of Dy , X t , X^): 

32 1 

-g 2 3 M 2 + -g 2 D Y , (4.36a) 

6g 2 M 2 - ys§Mf + \g\D Y . (4.36b) 

Note that in the RG equations of all soft masses, dependence on the gaugino mass 
Mi occurs only as terms proportional to g\M."^. Hence we can eliminate the M2 
dependence by taking suitable linear combinations of MSSM parameters. For exam- 
ple, the RG equation for the quantity m 2 ~ — m 2 ~ (which occurs in Db 13 ) does not 

t^l W 3 

depend on M2 any more, so its experimental uncertainty will spread less under RG 
flow. However, in this process we have thrown away information about the value of 

m ~ + m 2 ~ . Thus we have to reduce the number of independent quantities to reduce 

Qi Q'a 

the spread of uncertainties under RG flow. 

This may become a limitation of the RGI method in the following sense. A minimal 
model such as mSUGRA, with only three parameters that govern the soft masses plus 
gauge couplings at the GUT scale (mo, M1/2, g = g a (^GUT)), allows us to construct 
sum rules because we have more RGIs than mSUGRA has parameters. However, if 
we have a not-so-minimal model with (say) 15 parameters that determine the high- 
scale spectrum, we do not have enough RGIs to make any sum rules. 7 Hence, despite 
the simplicity of the method, we are still limited to models with few parameters. 

• The applicability of the RGI method to the study of supersymmetry breaking depends 
crucially on the assumption that the MSSM renormalisation group equations are valid 
all the way up to the scale of supersymmetry breaking. But suppose that in Nature 
a new field <I> (or possibly more than one) enters the theory at a high scale 
that is not the scale of supersymmetry breaking; instead supersymmetry is broken 
at an even higher scale hsvsy- Then at /i$ the physical RG trajectories of the 
running parameters will be deflected from their MSSM trajectories. Thus we might 
mistakingly see gaugino-mass unification where it is absent, or vice versa. Hence, if 
we want to study supersymmetry breaking directly from RGIs, we have to assume 
that new physics, if present, does not alter the one-loop RG equations for the MSSM 
up to the scale of supersymmetry breaking. 

• In order to make conclusive statements based on sum rules, the values of all RGIs 
should be reconstructed. To achieve that, all soft masses and gauge couplings need to 
be known at one energy scale. This may prove difficult in practice. First of all, due to 

7 That is, unless a subset of n RGIs accidentally depends on less than n model parameters. 
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mixing effects the gauge eigenstates do not always correspond to the mass eigenstates. 
Reconstructing the soft masses from measured pole masses will introduce additional 
uncertainties. Furthermore, determining all soft masses and gauge couplings is one 
thing, but determining all of them at the same energy scale may prove challenging. 
Note however that the bottom-up method also suffers from these complications. 

• If we find some sort of unification to be consistent with experimental data, the RGI 
method does not provide us with any information on the value of the new-physics 
threshold. In order to determine the value of the scale of new physics, we need to 
resort to other methods such as the bottom-up method. 

5 Conclusions and Outlook 

We advocate to employ sum rules in terms of renormalisation group invariants as a sim- 
ple yet powerful method to probe high-scale physics in lower-energy experiments. This 
method has been worked out in detail for the study of supersymmetry-breaking mecha- 
nisms in the context of the Minimal Supersymmetric Standard Model. It has been argued 
that important clues about the supersymmetry-breaking mechanism are to be found in 
patterns between the high-scale soft-supersymmetry-breaking parameters. The renormali- 
sation group is the prime tool to extract such information on the high-scale spectrum from 
lower-energy data. Several methods have been discussed to do this and a new strategy 
has been proposed to make effective use of renormalisation group invariants. Assuming 
that the Minimal Supersymmetric Standard Model is an appropriate effective field theory 
beyond the Standard Model, a model-independent set of renormalisation-group-invariant 
sum rules has been constructed that test properties that are common in supersymmetry- 
breaking models, such as unifications and flavour-universality. If a certain property is 
realised in Nature, all corresponding sum rules must be satisfied. Since none of these sum 
rules refer to any parameters that are specific to some supersymmetry-breaking mechanism, 
they are useful regardless of the way supersymmetry has been broken in Nature. 

In addition, sum rules that are tailor-made for testing specific supersymmetry-breaking 
mechanisms have been considered. Their primary use was to determine the effectiveness 
and ambiguities associated with the model-independent sum rules. It was found that some 
sum rules do not provide unambiguous checks by themselves; however, in almost all cases 
other sum rules lift the ambiguity. Hence, for the currently known supersymmetry-breaking 
mechanisms the proposed model-independent sum rules are surprisingly effective. In the 
exceptional case when they are not, one may have to resort to other methods to resolve the 
ambiguity, such as a bottom-up analysis. It is possible that new breaking mechanisms will 
be proposed in the future, and that their corresponding sum rules introduce new ambigui- 
ties. Therefore, it should be kept in mind that the main strength of invariant sum rules is 
their falsifying power. If we are able to determine all soft masses and gauge couplings, the 
compatibility of the sum rules with experimental data will put severe constraints on any 
realistic model of supersymmetry breaking . 

It is possible that the next effective field theory beyond the Standard Model is not the 
Minimal Supersymmetric Standard Model. It may as well be a non-minimal supersymmet- 
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ric extension of the Standard Model, or even a non-super symmetric theory. Nevertheless, 
the proposed scheme for probing high-scale properties of running parameters may be ap- 
plied just the same. In order to perform an analogous study, one needs to determine 
the particle content, interactions and /3- functions of the appropriate effective field theory. 
Then one should determine all independent renormalisation group invariants for this effec- 
tive field theory and construct sum rules in a way similar to what has been worked out in 
this study. However, there is a large amount of structure in the /^-functions of the Minimal 
Supersymmetric Standard Model in view of the limited number of combinations in which 
running parameters appear in them. So, an interesting topic is to determine the form of the 
(one-loop) /3-functions for a theory that is more general than the Minimal Supersymmetric 
Standard Model, and to see what renormalisation group invariants can be found for such 
a general theory. We leave these issues to future work. 

Apart from establishing the validity of the Minimal Supersymmetric Standard Model, 
the main obstacle to using our invariant sum rules is the necessity of knowing all soft 
masses and gauge couplings at one scale. Therefore, an important topic for future study 
will be to determine how well this can be done and how the sum rules can be exploited 
as a fast means of identifying theoretically interesting regions in as yet unconstrained 
parameter space. Another important issue is to find out how the soft mass parameters can 
be reconstructed from the mass eigenstates of the sparticles. 
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A One-loop RG equations for the MSSM 



In this appendix we give the renormalisation group equations of the MSSM that have been 
used in this study. They have been taken from [17] and are one-loop equations that have 
been simplified by the assumptions for the pMSSM (see section 3.3). For general two-loop 
RG equations, see e.g. [38]. It is convenient to use the /^-functions, which differ from their 
corresponding RG equations by a constant: 

m = i 6vr 2^. (A.!) 

Here p is a running parameter and t = log (fi/fio), where /i is the renormalisation scale and 
/xo an (arbitrary) energy scale that makes the argument of the logarithm dimensionless. 
Under the approximations of the pMSSM we are left with the following running parameters: 

g a (a = 1,2,3) Gauge couplings 

M a (a = 1,2,3) Soft-supersymmetry-breaking gaugino masses 

m ~ , m| , ml , m~ , m| Soft-supersymmetry-breaking sfermion masses 

m 2 Hu , vr?n d Soft-supersymmetry-breaking Higgs mass parameters 

ytiVbiVT Yukawa couplings for the third-generation (s)fermions 

Af,Ab,A T Soft-supersymmetry-breaking trilinear couplings for the 

third-generation sfermions 

[i Supersymmetry-respecting Higgs mixing parameter 

B Soft-supersymmetry-breaking Higgs mixing parameter 

Here we use the soft Higgs mixing parameter B = b/[i rather than b because its /3- function is 
simpler. For the sfermion masses, we denote the first and third generation with a subscript 
1 and 3 respectively. 

The following notation is used for parameters that enter the RG equations through 
common combinations of Dynkin indices and quadratic Casimir invariants: 

b a =(j,l,-3\ for a = 1,2, 3. (A.2) 

It is also convenient to define the following combination of running parameters, which 
appears in the RG equations of the sfermion masses: 

Dy = Tr (7m 2 ) = ( m <5 — m | — m Z m l) — m #d' (A-3) 

gen 

Here the trace runs over all chiral multiplets and the sum runs over the three sfermion 
generations. Note that Dy is often called S in the literature. Furthermore, we define the 
useful combinations: 

X t = 2\y t \ 2 (m 2 Hu + m| s + m| 3 + \A t \ 2 ) , (A.4a) 

X b = 2\y b \ 2 (m 2 Hd + m| a + ml + |A b | 2 ) , (A.4b) 

X T = 2\y T \ 2 (m 2 Hd + m| a + m| + |A r | 2 ) . (A.4c) 
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Then the resulting /3-functions for the MSSM are: 
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(A.5a) 
(A.5b) 

(A.5c) 
(A.5d) 
(A.5e) 
(A.5f) 
(A.5g) 
(A.5h) 
(A.5i) 
(A.5j) 
(A.5k) 
(A.51) 
(A.5m) 
(A.5n) 
(A.5o) 
(A.5p) 
(A.5q) 
(A.5r) 
(A.5s) 

(A.5t) 
(A.5u) 
(A.5v) 
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B Deriving the one-loop RGIs for the MSSM 



In this appendix we will derive a maximal set of independent RGIs for the MSSM. First 
we will determine invariants that contain the running parameters fi and B = b/fi. We 
will see that there is only one independent RGI for each of them, making them useless for 
our study. Then we will argue that we are restricted to RGIs containing only soft masses 
and/or gauge couplings. We will derive all of them systematically; our approach will be 
globally the same as in [8], but using different arguments to show that we do indeed find 
all RGIs. 

Let us consider the parameter fx. The only /3-function containing fi is that of fi itself. 
Note that we can write j3{fi) more conveniently as: 



/3(log/i) 



3M 2 + 3|y fe | 2 + |y T | 2 



-9\ 



3gh 



The only other /3- functions containing terms linear in \yt\ 2 , |j/b| 2 , \ Ut 
logarithms of the Yukawa couplings: 



Pi}ogy t ) 
£(log Vb) 
/3(logy r ) 



Q\yt\ 2 + \y b \ 2 



13 2 16 2 

— g 1 - 3# 2 - y£ 3 , 

I I o r>\ \1 I 1 9 ^2 o 2 2 

\yt\ +6|y b | + \y T \ - — g 1 - 3g 2 - y# 3 , 

3\y b \ 2 + A\y T \ 2 - Igj - 3g 2 2 . 
5 



(B.l) 

are those of the 

(B.2a) 
(B.2b) 
(B.2c) 



The terms in the /3-functions proportional to g 2 can be eliminated by taking linear com- 
binations with logarithms of gauge couplings, of which we can rewrite the /3-functions 
as: 

/3(log 5 a) = b a g 2 a (a = 1,2,3). (B.3) 

Hence, fi can only appear in an RGI through a linear combination of log/i, logyt, logy;,, 
logy T , logyi, log<?2 and log 53. 8 We have seven /3-functions with six different terms to 
eliminate (namely terms linear in \yt\ 2 , \yb\ 2 , \yr\ 2 , g 2 , g 2 or g 2 ), so we can make one RG 
invariant linear combination of them. Using elementary linear algebra we find that the 
linear combination 
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- 77 log Vt - 77 lO; 
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61 
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- 77 l°g Vr + log fi - 77 • 77 log 51 

61 33 
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9 , 1 256, 

+ 77 log 92 + — ■ — log 53 
61 61 3 

9 256/3 1 V61 
92 93 



yfyfy^g[^\ 



(B.4) 



We could also include logarithms of gaugino masses in these linear combinations, since their /3-functions 
are also proportional to g\. However, in a moment we will construct RGIs from the gauge couplings and 
gaugino masses only. Any RGI that contains both \i and the gaugino masses will be a function of those 
RGIs and the one we are constructing now. 
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has a vanishing /3-function. Thus we can choose the only independent RGI containing \i to 
be 



q 256/3 
92 9 3 



yfyfy^g[^\ 



1/61 

(B.5) 



using the notation of [7] . To summarise, we have found a set of independent RGIs containing 
\i (in this case only one) by considering what terms in the MSSM /3-functions could cancel 
each other. This will be our general strategy for finding all RGIs of the MSSM, because the 
running parameters only enter the /3-functions in a very limited number of combinations 
(e.g. the soft scalar masses only appear in the linear combinations Dy, Xt, -X& and X r ). 

Now we turn to the parameter B. It does not appear in any of the MSSM /3-functions 
itself. Its /3-function contains only terms linear in ^4^ | j/t | 2 , ^4fo| j/t,| 2 , A r |y T | 2 , g\Mi, g\Mi 
and g\M%. The /3-functions of At, A^, A T , Mi, Mi and M3 also contain only these terms, 
so B should always appear in RGIs in a linear combination of these parameters. This gives 
us seven /3-functions with six different terms to eliminate, so again we can make one RG- 
invariant linear combination. Using elementary linear algebra this combination is found to 

4 = B - >V- J* + JLJ,,. (B.6) 
61 61 61 183 61 2013 v ' 

Indeed we have found only one independent RGI containing \x and one containing B. 
As was argued in section 3.10, RGIs are only useful as long as their constituent running 
parameters also appear in other RGIs. This is not the case for I2 and I4, so we are restricted 
to RGIs that contain neither /i nor B. But in the above procedure, we needed their /3- 
functions to eliminate the \yi\ 2 and j4j|yj| 2 dependence respectively from the /3-function of 
the RGI under construction. If we wish to construct RGIs containing the Yukawa couplings 
without using [i, we have to eliminate three different |yj| 2 terms using three /3-functions, 
so we cannot make any RG-invariant combinations. Similarly, we cannot make any RGIs 
containing the soft trilinear couplings without using B, because we have to eliminate three 
different Aj|?/j| 2 terms using three /3-functions. 

Thus, if we want to construct RGIs without using fi and B, we cannot use the Yukawa 
and soft trilinear couplings either: we do not have enough equations to eliminate all terms 
from the /3-function of the RGI under construction. Therefore, from now on we will only 
consider RGIs that are functions of soft masses (fifteen parameters) and / or gauge couplings 
(three parameters). 

Let us begin with RGIs constructed from the gauge couplings only. First we rewrite 
their /3-functions into a more convenient form: 

66 

P(9i 2 ) = -y , (B.7a) 
(3(g 2 2 ) = -2, (B.7b) 
P{gf) = 6. (B.7c) 

This gives us three equations to eliminate a single term (namely a constant), hence we can 
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make two independent RGIs out of them. In accordance with [8], we choose them to be 



I 92 ^g^-fg^ 2 , (B.8) 
I g3 =9i 2 + ^9z 2 . (B.9) 

Now we turn to the gaugino masses. First we rewrite their /3-functions as follows: 

66 

/3(logM l ) = -gl (B.lOa) 
5 

/3(logM 2 ) = 2g 2 2 , (B.lOb) 
0OogM 3 ) = -6g§. (B.lOc) 

Together with (B.3) this gives six equations with three different terms (namely those pro- 
portional to g 2 ) to eliminate. Hence, we get three new RGIs by taking linear combinations 
of logM a and \ogg a : 

(B.lla) 
(B.llb) 
(B.llc) 

93 

Thus we can choose the three independent RGIs to be: 

I Bl = ^, (B.12) 
9\ 

Ib 2 = — 5") (B.13) 

" si 

IBs = ^- (B.14) 

93 

Now let us consider RGIs constructed solely from the twelve soft scalar masses. First we 
eliminate the Yukawa terms Xt, Xj,, X T and the gaugino-mass terms g 2 M 2 , g^M^, g 2 M 2 
from the /?- function. Since we have to eliminate six terms using twelve equations, we can 
make six independent linear combinations of the soft scalar masses that have a /3-function 
proportional to g\Dy. Then we can make linear combinations of these quantities such that 
five of them have a vanishing /3-function and the sixth quantity still runs with g\Dy . In 
accordance with [8], we choose the five RGIs to be: 9 

Dn,,, = 2 ( m 2 ~ — m 2 ~ ) — m~ + m~ — m~ + m~ , (B.15) 
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= /3(log 


M 3 
1 2 








9s 



9 The notation used for the RGIs may look odd here. In [8] they are related to symmetries of the MSSM 
Lagrangian. In this context, the D-term Di of a charge Qi is defined as Di = Tr(Qim 2 ), with the trace 
running over all chiral multiplets. Then one should interpret Db 13 as Db 1 — Db 3 , where the subscripts 1 
and 3 mean that the trace is restricted to sfermions of the first and third generation respectively. See [8] 
for an explanation of the nomenclature for the remaining RGIs. 
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D Ll3 


= 2 (ml 




= 3 f 3m ; 


D Y 1SH 





-m~ ) -ml + ml , (B.16) 

l LzJ ei e 3 ' V > 

I -2(m% -ml) -ml) -ml , (B.17) 
m ~ — 2m~ + ml — m~ + ml 

Qi «i d! L! ei 

10 / 2 r,2,2 2,2,2 2 \ / D10 \ 

m~ — 2m~ + m^ — m r + m~ imri — m ffj , B.lo 

D z = 3 (m| - m| ) + 2 ^m| g - mfj . (B.19) 

The sixth quantity, which runs with gfDy, can be chosen to be Dy itself, because 

66 

(3{D Y ) = -gfDy. (B.20) 
5 

Note that log Dy runs with gf, so using (B.3) we find 

/3(lo gj Dy-21og 5l ) = /3(log^) = 0. (B.21) 

9i 

This gives us another independent RGI: 

fr. = = 4 ( ™k - + £ (m| - 2m\ + m| - m\ + mf) ) . (B.22) 

"1 "1 \ gen 



Finally, we look for RGIs constructed from both scalar masses and gaugino masses. Note 
that the gaugino-mass /3- functions can be rewritten as: 

(3(M 2 a ) = ^ a g 2 a Ml (a = 1, 2, 3). (B.23) 

Combining the gaugino masses and scalar masses, we have fifteen /3-functions with seven 
terms to eliminate, so we can construct eight RGIs by taking linear combinations of the 
squared gaugino masses and the scalar masses. Five of them can be made from the scalar 
masses alone, so there must be three new RGIs. In accordance with [8], we take them to 
be 

I Ml = Mf - | (m| - ml - m|J , (B.24) 
h h =M% + ± (9 (m| - m|J + lBm 2 ^ - m|) , (B.25) 
I M3 ee Ml - A (5ml + ml - m|) . (B.26) 
These complete the list of independent one-loop RGIs for the MSSM. 
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